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We report on the residence times of capillary waves above a given 
height h and on the typical waiting time in between such fluctuations. The 
measurements were made on phase separated colloid-polymer systems by 
laser scanning confocal microscopy. Due to the Brownian character of the 
process, the stochastics vary with the chosen measurement interval At. 
In experiments, the discrete scanning times are a practical cutoff and we 
are able to measure the waiting time as a function of this cutoff. The 
measurement interval dependence of the observed waiting and residence 
times turns out to be solely determined by the time dependent height- 
height correlation function g{t). We find excellent agreement with the 
theory presented here along with the experiments. 

1 Introduction 

The often counterintuitive field of stochastic number fluctuations has boggled the 
scientist's mind for ages. Svedberg's measurement of the diffusion constant [1], based 
on the residence times on a fluctuating number of Brownian particles in a fixed re- 
gion, inspired von Smoluchowski to develop a theoretical framework for it [21 E] • His 
ideas were for example exploited to derive the mobility of spermatozoids |1] and white 
blood cells [5]. The original experiment was performed in even more detail by Bren- 
ner, Weiss and Nossal |6]. Based on Einstein's [7] and Perrin's seminal papers on 
Brownian motion of particles [8], von Smoluchowski |9] was the first to predict the 
Brownian height fluctuations of the interface. They were first theoretically treated 
by Mandelstam [TU] and have become an important component of modern theories 
of interfaces pTl [T2| 113] . Capillary waves were initially accessed experimentally by 
light- [H] and x-ray scattering [16] . On a microscopic level capillary waves were 
studied in computer simulations of molecular systems [17|, before recent investiga- 
tions by Aarts and coworkers on colloid-polymer mixtures, [T9| [20] added another 
dimension to studies on capillary waves by using confocal microscopy. In these ex- 
periments, the interfacial tension 7 is lowered to the nN/m range. As a consequence 
the characteristic length and time scale of the fluctuations are such that they can be 
visualized by microscopy. Microscopy furthermore enables the investigation of local 
phenomena such as the role of capillary waves in the rupture [21] and coalescence [22] 




Figure 1: (A) Variation of the interfacial height h in time with characteristic residence 
times T and waiting times 0. (B) in practice, T and are determined by discrete 
intervals, resulting in observed residence and waiting times and r^. These are 
determined by weighted averages, of the stochastics of intervals with length At. 



problems. At a fixed location on the interface, the height is continuously fluctuating 
in time. Following Becker [23j and as sketched out in Figjl]\ the waiting time G is 
the average time spent in between fluctuations above a height /i, and the residence 
time T is the average time spent above a height h: 



Q{h) 



dtp-{h,t)t, T{h) = dtp+{h,t)t 



(1) 



where p~{h,t) and p~^{h,t) are the probabilities on intervals of length t below and 
above the height h respectively. Such local measurements are not possible by scatter- 
ing methods, since they require knowledge of continuous local stochastics which seem, 
however, easily accessible by microscopy. Pioneering work by Aarts and Lekkerkerker 
[21] resulted in scaling relations, but did not include a full quantitative description of 
the process. The crux is that however fast microscopy may be, measurements always 
have to be taken at discrete time intervals At, as sketched in Fig. [Tj3, resulting in 
the observed waiting and residence times r+ and r~: 



(2) 



Here, p~ and p^ are respectively the (normalized) probabilities on intervals of snap- 
shots of length n, below and above height h. Note that and r~ are in units At and 
are therefore both dimensionless. Switching to discrete time intervals is not entirely 
trivial. Due to the Brownian character of the process, the discretisation of [T] leads 
to statistics that depend on the chosen interval At. However, the necessity to be 
discrete saves rather than spoils the day, as it enables us to overcome the divergencies 
for continuous distributions. 
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Figure 2: (A) Phase diagram of the system studied. The dilution hne (hne 11) studied 
here is marked. All shown state points phase separate. The indicated binodal is a 
guide to the eye: the theoretical binodal occurs at much lower concentrations. The 
marked critical point is an estimate based on the ratios of the volumes of the "liquid" 
and "gas" phases. (B) Confocal Images of Statepoint II and IV. 
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2 Experimental Section 



Fluorescently labeled polymethyl-metacrylate particles were prepared using the Bosma 
method [25] , slightly modified by using decalin (Merck, for synthesis) as the reaction 
solvent [26j. The particle polydispersity is around 10% from scanning electron mi- 
croscopy, and the dynamic light scattering particle radius is 62 nm. Polystyrene (2 
10^ kg mol~^) was added as depletant polymer, with an estimated radius of gyration 
Rg of 42 nm [27]. At sufficiently high colloid and polymer volume fractions, respec- 
tively (pc = 7r/6(j^nc and 0p = i/SnRgnp (with ric and Up the number densities of 
colloids and polymers respectively), this system phase separates into a colloid-rich 
(colloidal liquid) and a polymer rich (colloidal gas) phase. By diluting several phase 
separating samples with its solvent decalin, the phase diagram presented in Fig. |2]A. 
was constructed. The shown binodal is a guide to the eye (the theoretical binodal 
appears at much lower volume fractions) and the critical point is an estimate based 
on the ratios of the volumes of the phases. 

A Nikon Eclipse E400 laser scanning confocal microscope equipped with a Nikon 
CI scanhead was placed horizontally to study the colloid polymer mixture [26]. The 
microscope was furthermore equipped with a 405 nm laser and a Nikon 60X CFI 
Plan Apochromat (NA 1.4) Lens. The sample container is a small glass vial, part 
of which is removed and replaced by a thin (0.17 mm) glass wall. Series of 10 000 
snapshots of the interface of 640 x 64 and 640 x 80 pixels (100 x 10.0 and 12.5 fim) 
were taken at constant intervals ti of 0.45 s and 0.50 s of statepoints II and IV along 
the marked dilution line in the phase diagram, shown in Fig. [2j A single scan takes 
approximately 0.25 s to complete (the exact scan time does vary by a few percent in 
time). Typical snapshots are shown for a number of state points in Fig. [2]d. The low 
excitation wavelength results in resolution of ~ 160 nm. The particles are ~ 138 nm 
in diameter, hence a pixel roughly corresponds to a particle. Note that the resolution 
of the measured heights is significantly higher than this: the vertical location of the 
interface h{x) is determined for each column of pixels in a frame by fitting the pixel 
value I{z), which is proportional to the local colloid concentration, to a van der Waals 
profile I{z) = a + b tanh([z — h{x)]/c), as in [19j. The exact resolution depends on 
the contrast between the phases, which depends on the distance to the critical point. 
Further away from the binodal the height fiuctuations are too small, closer to the 
binodal the contrast between the phases vanishes and the capillary waves start to 
show overhang effects. 



3 Observations and Results 

The wave heights h at an arbitrary point are distributed according to the Gaussian 
It can be shown that {h'^) satisfies the expression 



T2 pK,nax U 7^ J- 

^ / /IL |2\; 77. '^Bi^ i.'^max + S 
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Figure 3: Interfacial properties determined by dynamic correlation functions and 
height distributions for state points II and IV. (A) height distributions for several 
state points. (B) The dynamic correlation functions obtained for statepoints II and 
IV, normalized by (/i^). 

Here ^ = {'j/ ApgY^'^ is the the capillary length, with g the gravitational acceleration, 
Ap the density difference between the two phases, the Boltzmann constant, T 
the Kelvin temperature and L the system size. Furthermore, hk is the amplitude of 
mode k, k^ax ~ 2Tr/d is a cutoff related to the typical interparticle distance d and 
kmin = 271 /L is related to the system size L. In equation Q we have used that the 
mean square average of hk is given by 



The height distributions are shown for statepoints II and IV in Fig. [sjA., with (/i^)^/^ = 
0.219 and 0.336/im for statepoints II an IV respectively. These will be used as a unit 
for the heights. 

Next we consider the height correlations in time {h{t)h{t')) at a fixed position, 
which we denote by {h'^)g{t — t'). It is calculated as: [28l 129] 

{h^)g{t-t') = ^ (|/ifcHe""^l*^*'lA;rfA;. (6) 



2rr Jt 

The measured and fitted (setting k 

max — CK) and fcmm — 0) dynamic correlation 
functions of statepoints II and IV are shown in Fig. [3^. In colloidal systems capillary 
waves are in the overdamped regime [281 [H] with a decay rate 

u^k = ^^{k^+{kO-')- (7) 

Here, = ^/uc is the capillary time and = j/f] is the capillary velocity with rj 
the combined viscosities of the phases. From the fits, the interfacial tensions and the 
capillary times can be extracted, which gives tc = 12s and 7 = 69nN/m for statepoint 
II and tc = 22s and 7 = 25nN/m for statepoint IV. 
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Figure 4: (A) The distribution of observed residence times r"*" for statepoint IV. 
Pn{h) is shown for h = -1 (filled squares), (circles), and 1 (semi-filled pentagons) 
X The inset shows the short time behaviour for /i = in detail. The colours 

correspond to time intervals A t = 1 (black squares), 2 (red circles) and 4 tj 
(green triangles). (B) Observed waiting times (open symbols) and residence times 
(filled symbols), observed for statepoint II, at time intervals of 1 ti (black squares), 2 
ti (red circles) and 4 ti (green triangles). (C) as in panel (B), but now for statepoint 
IV. 

The residence and waiting times are calculated through the experimentally ob- 
tained functions p^{h). In order to check the possible variation of waiting and res- 
idence times, we calculate the distributions for intervals of 1,2 and 4 ti. They are 
shown for statepoint IV at heights h = —1,0 and l(/i^)^/^ in Fig. |4]A.. Note that 
the X-axis has intervals as units, not time. The distributions are clearly complicated: 
the shortest interval of 1 time interval becomes more dominant as At decreases as 
shown in the inset of Fig. On the other hand, h = — (/i^)^ has a surprisingly 
large interval of ~ 2-50 interval durations, which occur most frequently for At = 4tj. 
The longest intervals decay exponentially and as expected decrease most rapidly in 
quantity for the largest time intervals. As a result, the waiting and residence times, 
shown for statepoints II and IV as a function of height in Fig. |4]B and Fig. |4p, are 
confusing at first sight: if we increase the time in between observations At to 2 or 4 
ti, the observed residence times are significantly larger! 
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4 Interpretation and Discussion 



In order to understand these observations, we need some theoretical framework. In 
the spirit of von Smoluchowski we start with simple counting arguments, [2l E], in 
order to identify the essential object to be calculated. Consider a long measurement 
of snapshots. We take pairs of consecutive heights and divide them into 4 sets: 
N~^'^{h) pairs where on both sides of the interval the interface is above h, N^~{h) 
pairs where the earlier value is above h and the later value below h, N~~^{h) pairs 
where it is the other way around and finally (h) pairs where the interface is on 
both sides of the interval below h. The first observation is that 

N+-{h) = N-+{h) = M{h), (8) 

since each interval in which the interface crosses the level h from above is followed by 
the next crossing from below. These numbers are also equal to the number of M{h) 
of "hills", which in turn is the same as the number of "valleys". Now M{h)T~^{h) is 
the total length of the hills (measured in units At). Thus we have the relations 

M{h)T+{h) = N++{h) + M{h), M{h)T~{h) = N''{h) + M{h). (9) 

We have to add M[h) on the right hand side, since the number of points in a hill 
(valley) is one more than the number of intervals inside a hill (valley). Adding the 
two relations gives 

-"<'" = sw"Rb- 

r{h) is the probability to find an interval where the interface crosses the level h from 
above to below h. 

The second observation is that the right hand side of ^ gives the probability to 
find a point above (below) h\ 

N++{h) + M{h) , r ,,,r, N-~{h) + M{h) 



N 



oo 

fill 



which can be calculated from the equilibrium distribution and which states that the 
fraction of heights above or below h is equivalent to the fraction of time spent above 
or below height h. Thus we find the expressions for T^{h) 



showing that r{h) is the basic quantity to be calculated. Equation ( 12 ), independently 
of r{h), implies: 

T-ih) q-ihy ^ 

which is experimentally verifiable. In Fig. |5]A we show that both sides of 13 fall on 



a single mastercurve over nearly the full spectrum of heights for both statepoints. 
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For r{h) we need the joint probability to find the interface on level h' at t = and 
on h" at time t 

{S{h{0, 0) - h')6{h{0, t) - h")) = P,^{h')Gc{h', 0\h", t). (14) 

It is the product of the height distribution function and the conditional probability, 
that starting from h' a.t t = one arrives at h" a time t later. Gc follows from the 
Langevin equation for the interface modes k, using a fluctuation force with a white 
noise spectrum [30]. The result is 



We point out that Gc contains only the correlation function g{t) and not explicitely 



the decay rate u^- A check on (15) is that the correlation function {h{0)h{t)) as 



calculated from (14) and (15) gives the result {h'^)g(t). Fig. [Sp shows that g(t) is 
not exponentially decaying. Thus the spatial process is non-Markovian [30] . 
According to 14 r{h) is given by the integral 

POO ph 

r{h,At)= dh'Pegih') dh"Gcih',0\h",At), (16) 

J h J — oo 

where the first integral selects the points above h with the height distribution and 
the second integral selects the points below h with the conditional probability given 



by (15). We have made the dependence on the interval length At explicit because 
it is an important issue. At enters in the calculation of r{h,At) through the value 
g{At). We may therefore consider r also as a function of h and g = g{At), using the 
experimental values of g{t) as input. Of course r{h,g) is experimentally measurable 



using its definition (10). The theoretical and experimental r(h,g) are plotted in Fig. 
|5p as function of g for a few values of h for both statepoints, resulting in a single curve 
for each h/ {h^)\- The agreement between theory and experiment is remarkable. The 
observed residence and waiting times for At = 1,2 and Ati for statepoints II and IV, 
are shown in Fig. [5p and[5p in terms of units At. Again we find excellent agreement 
with the theory, using the measured values for the g{t) function as theoretical input 
for the determination of r. Note that the y-axis has units "At" and is dimensionless. 
In terms of time, the residence time at height h increases with At. At sufficiently large 
the theoretical times are always larger than the experimental values, see Fig. [5p 
and [Sp. This is at least partly due to the limited amount of time points considered, 
which tends to exclude the tail of the distributions Pn{h). 

We now ask ourselves the question whether we can take the limit of At — >■ and 
arrive at a continuum description. If r{h, At) were to shrink proportional to At, the 
T^{h) would increase inversely proportional to At. From a continuum description one 
would expect the relations 

T = lim r+(/i)At, = lim r"(/i)At (17) 

At-»0 At^O 

and T and B would have a finite limit. However, analyzing the small At behavior of 



the expression (16), which amounts to the limiting behavior for g{At) — > 1, we find 



r{h,g ^ 1) ~ — — ^exp 



h 



2 
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Figure 5: Waiting and residence times for statepoints II and IV. (A) The relation 
q^{h)/q^{h) = T^{h)/T^{h) clearly holds over many orders of magnitude, both for 
statepoint II (dashed lines) and IV (non-dashed lines). (B) the probabihty r{h, g) as a 
function of g for \h\ = 0, 1 and 2 (/i^)^/^. The lines are the theoretical expressions, the 
points are obtained through the definition of r from the experiment. (C) Theoretical 
and observed time intervals for the waiting and residence times of statepoint II, with 
At = 1 (black squares), 2 (red circles) and 4 (green triangles) tj. Open symbols 
are waiting times, closed symbols are residence times. The lines are the theoretical 
curves. (D) Theoretical and observed time intervals for statepoint IV, with At = 1,2 
and 4 tj. Symbols are identical to those in panel (C). 
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The square root in ( 18 ) is a reflection of the Brownian character of the fluctuations. 
If the initial decay of g{t) is hnear, r{h,g) vanishes as (At)^/^ and consequently T 
and vanish. 

With 12 and the computation of r{h,g) we have determined the residence and 
waiting times without evaluating the distributions Pn{h). The latter quantities can 
also be calculated using multiple correlation functions, which are generalizations of 



14} Their evaluation is increasingly involved, due to the non-Markovian character of 
the spatial process. Experimentally, the exhaustive amount of multiple correlation 
functions can be determined given sufficient time points, analogous to the results in 
Fig. |5j3. We leave this aspect for a further study. 



5 Conclusion 

We have presented confocal microscopy experiments along with theory for the micro- 
scopic waiting and residence times of heights h of the capillary waves of the fluid-fluid 
interface of a phase separated colloid-polymer mixture. Due to the Brownian charac- 
ter of the process, these times depend on the experimental measurement interval At. 
The results from this discrete time sampling are predictable in terms of the decay of 
the height-height correlation function g(t), which in principle includes the relevant 
system variables interfacial tension and capillary time. We flnd excellent agreement 
between experiments and theory. 

Acknowledgements The authors are indebted for valuable discussions to 
H.W.J. Blote and W.Th.F. den Hollander. We additionally thank C. Vonk for assis- 
tance with particle synthesis, S. Sacanna for help with the electron microscopy and B. 
Kuipers for aid with the dynamic light scattering experiments. The work of VWAdV 
is part of the research program of the 'Stichting voor Fundamenteel Onderzoek der 
Materie (FOM)', which is flnancially supported by the 'Nederlandse Organisatie voor 
Wetenschappelijk Onderzoek (NWO)'. Support of VWAdV by the DFG through the 
SFB TR6 is acknowledged. 

References 

[1] T. Svedberg. Z. f. Phys. Chem., 77:147, 1911. 

[2] M. Von Smoluchowski. Wien. Ber., 123:2381, 1914. 

[3] M. Von Smoluchowski. W%en. Ber., 124:339, 1915. 

[4] L. Rothschild. J. Exp. BioL, 30:178, 1953. 

[5] R. J. Nossal and Y. T. Chang. J. Mechanochem. Cell Motii, 3, 1976. 

[6] S. L. Brenner, R. J. Nossal, and G. H. Weiss. J. Stat. Phys., 18:1-18, 1978. 

[7] A. Einstein. Ann. Phys., 17:549, 1905. 



10 



[8] J. Perrin. Ann. de Chim. Phys., 18:5, 1908. 

[9] M. V. von Smoluchowski. Ann. Phys., 25:205, 1908. 

[10] L. Mandelstam. Ann. Phys., 41:609-624, 1913. 

[11] F. P. Buff, R. A. Lovett, and F. H. Stillinger. Phys. Rev. Lett, 15:621, 1965. 

[12] K. R. Mecke and S. Dietricli. Phys. Rev. E, 59:6766, 1999. 

[13] A. Milcfiev and K. Binder. Europhys. Lett, 59:81, 2002. 

[14] A. Vrij. Adv. Coll. Interf. Sci, 2:39-64, 1968. 

[15] C. Fradin, A. Braslau, D. Luzet, D. Smiigies, M. Alba, N. Boudet, K. Mecke, 
and J. Daillant. Nature, 403:871, 2000. 

[16] M. K. Sanyal, S. K. Sinlia, K. G. Huang, and B. M. Ocko. Phys. Rev. Lett, 
66:628, 1991. 

[17] J.H. Sikkenk, J.M.J, van Leeuwen, E.O. Vossnack, and A.F. Bakker. Physica A, 
146A:622-633, 1987. 

[18] D. G. A. L. Aarts, M. Scfimidt, and H. N. W. Lekkerkerker. Science, 304:847, 
2004. 

[19] D. Derks, D. Aarts, D. Bonn, H. N. W. Lekkerkerker, and A. Imfiof. Physical 
Review Letters, 97, 2006. 

[20] C. P. Royall, D. G. A. L. Aarts, and H. Tanaka. Nature Physics, 2007. 

[21] Y. Hcnncquin, D. Aarts, J. H. van der Wiel, G. Wegdam, J. Eggers, H. N. W. 
Lekkerkerker, and D. Bonn. Physical Review Letters, 97, 2006. 

[22] D. Aarts. Soft Matter, 3:19-23, 2007. 

[23] R. Becker. Theorie der Wdrme. Springer- Verlag, Berlin, 1966. 

[24] D. G. A. L. Aarts and H. N. W. Lekkerkerker. J Fluid Mech., submitted. 

[25] G. Bosma, C. Pathmamanoharan, E. H. A. dc Hoog, W. K. Kegel, A. van 
Blaaderen, and H. N. W. Lekkerkerker. J. Colloid Interface Sci., 245:292, 2002. 

[26] D. G. A. L. Aarts and H. N. W. Lekkerkerker. J. Phys.: Condens. Matter, 
16:84231, 2004. 

[27] B. Vincent. Colloids Surf, 50:241, 1990. 

[28] U. S. Jeng, L. Esibov, L. Crow, and A. Steyerl. J. Phys.: Condens. Matter, 
10:4955, 1998. 



11 



[29] J. Meunier. Experimental studies of liquids at interfaces: classical methods 

for studies of surface tension, surface wave propagation and surface thermal 
modes for studies of viscoelasticity and bending elasticity. In J. Charvolin, 
J. F. Joanny, and J. Zinn- Justin, editors. Liquids and Interfaces, pages 327-369. 
North-Holland, New York, 1988. 

[30] N. G. van Kampen. Stochastic Processes in Physics and Chemistry. Elsevier, 
3rd edition edition, 2007. 



12 



